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ABSTRACT 

We study the biasing relation between dark-matter halos or galaxies and the under- 
lying mass distribution, using cosmological JV-body simulations in which galaxies are 
modelled via semi-analytic recipes. The nonlinear, stochastic biasing is quantified in 
terms of the mean biasing function and the scatter about it as a function of time, scale 
and object properties. The biasing of galaxies and halos shows a general similarity and 
a characteristic shape, with no galaxies in deep voids and a steep slope in moderately 
underdense regions. At ~ 8/i _1 Mpc, the nonlinearity is typically £ 10 percent and 
the stochasticity is a few tens of percent, corresponding to ~ 30 percent variations in 
the cosmological parameter (3 — fl 6 /b. Biasing depends weakly on halo mass, galaxy 
luminosity, and scale. The observed trend with luminosity is reproduced when dust 
extinction is included. The time evolution is rapid, with the mean biasing larger by a 
factor of a few at z ~ 3 compared to z — 0, and with a minimum for the nonlinear- 
ity and stochasticity at an intermediate redshift. Biasing today is a weak function of 
the cosmological model, reflecting the weak dependence on the power-spectrum shape, 
but the time evolution is more cosmology-dependent, relecting the effect of the growth 
rate. We provide predictions for the relative biasing of galaxies of different type and 
color, to be compared with upcoming large redshift surveys. Analytic models in which 
the number of objects is conserved underestimate the evolution of biasing, while mod- 
els that explicitly account for merging provide a good description of the biasing of 
halos and its evolution, suggesting that merging is a crucial element in the evolution 
of biasing. 



1 INTRODUCTION 



The standard picture of the growth of structure via gravita- 
tional instability within a dark-matter dominated universe 
has led to a successful and predictive theoretical framework. 
However, in order to make direct contact with most observa- 
tions, the relationship between galaxies and the underlying 
dark-matter distribution must be understood. This relation- 
s hip has come to be loosely referred to as galaxy "bias ing 



Pellegrini 1998; Tegmark & Bromley 1998), or the evolution 
of biasing in time ( jSteidcl ct al. 1998| ; Adelbcrgcr et al. 1998 ) 



(Davis {it al. 1985; Bardccn ct al. 1986~Dckcl & KcesTM 



To the extent that the cosmological parameters are known, 
studying galaxy biasing will allow us to better understand 
the process of galaxy formation. Conversely, better under- 
standing of biasing is a crucial component in interpreting 
the results of methods that attempt to use galaxies as trac- 
ers of the underlying density field and thus dete rmine the 
cosmological density parameter Q (see reviews by Btrauss & 
Willick 19951 iDekel 1994 |Dekel, Burstein, fc White 1997j, 



Given the complexity of the process of galaxy formation, 
it would be surprising if the galaxy distribution traced the 
mass distribution in a simple way. Various physical mecha- 
nisms have been proposed that could lead to galaxy biasing, 
but the features of the biasing process remain highly uncer- 
tain, including, for example, non-linearity, scale dependence, 
and stochasticity. Even the direction of the biasing, i.e. bias 
or anti-bias, is sometimes uncertain. 

Yet, there are many observational indications that 
galaxy biasing does exist and is non-trivial; for example, 
t he dependence of galaxy clustering on type or environment 
flPressler 198o|; [Hermit et al. 1996|; |Willmer, Da Costa, & 



Dekel 1999) 



and 

Modern iV-body techniques have provided the basis for 
a great deal of progress in understanding the clustering of 
dark-matt er halos relative t o the underlying mat ter density 
field (e.g. Gross et al. 199£, Jenkins et al. 1998). Analytic 



approximations based within the hierarch ical-merger frame- 
work have led to a dditional insights ( e.g. Mo & White 1996 
and extensions, cf. Catelan et al. 1995). However, these halos 



do not necessarily correspond to galaxies. One reason for this 
is the "over-merging" problem: dark matter halos that are 
incorporated within a larger collapsed structure may lose 
their identities due to the limited numerical resolution of 
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the simulation. In addition, the biasing properties of the ha- 
los include only the effects of clustering due to gravity and 
do not reflect the astrophysical processes involving gas and 
stars, which are presumably important in determining the 
properties of real galaxies. 

Recent theoretical studies have used a var iety of tech- 



niques to ad d ress these issues. For example, Cen fc Os- 
triker (|l998|), [Katz, Hernquist, fc Weinberg (1998|) |Blan- 



ton et |l (1998D and [Blanton ct al. (1999J used cosmologi- 
cal simulations with hydr odynamics and ph enomenological 
recipes for star formation. |Colin et al. (1998 ) used very high- 
resolution dissipationless N-body simulations, in which the 
over-merging problem is largely overcome and there should 
be a good one-to-o ne correspondence between dark matte r 
halo s and galaxies. Narayanan, Bcrlind, fc Weinberg (199£ 



and Mann, Peacock, fc Heavens (1998|) used various "toy 



models" to describe galaxy formation in large, low resolu 
tion dis sipationless N -bod y simulations 



Kauffmann, Nusser 



& Steii.metz (1997) and Benson et al. (199!:) used semi 



analytic techniques to assign galaxies to halos within dissi- 
pationless N-body simulations. 

All of these techniques have certain advantages and dis- 
advantages. The present study provides a useful complement 
to this previous work, and is unique in se veral respects. We 



modelgalaxies using a new technique (Kauffmann et al 



1998a|), in which iV-body simulations are combined with 
semi-analytic modelling of galaxy formation. Merging his- 
tories of dark matter halos (hereafter referred to simply 
as halos) are extracted from large dissipationless N-body 
simulations using the outputs at finely spaced time steps. 
Within the framework of these "merger trees" , gas dynam- 
ics, star formation, supernovae feedback, and galaxy merg- 
ing are modelled semi-analytically. The results are convolved 
with stellar population models. In this way, we obtain pre- 
dictions of observable galaxy properties such as magnitudes, 
colours, star formation rates, and morphologies. This allows 
us to study the dependence of galaxy biasing on these char- 
acteristics at different cosmological epochs and at different 
smoothing scales. Another novel aspect of this work is that 
we utilize the biasing formalism developed by Dekel & La- 
hav (l{f9g|), which allows us to separate and quantify the 
nonlinearity and the stochasticity of the biasing relation. 

The outline of the paper is as follows. In Section |^, we 
summarize the biasing formalism and analytic models for 
halo biasing. In Section ^, we briefly describe the simulations 
and the semi-analytic techniques used to model galaxy for- 
mation. In Section^, we discuss the relationship between the 
halos and galaxies in our simulations. In Section EL we study 
the dependence of the biasing relation on halo mass, galaxy 
luminosity and type, scale, and redshift. In Section m, we 
compare the results of the simulations with analytic mod- 
els of biasing. We summarize our results and conclude in 
Section |^. 



2 DESCRIPTIONS AND MODELS OF 
BIASING 

A biasing scheme relates the density fluctuation fields of 
galaxies and mass. For the mass, we define 8(x) = [p(x) — 
pb]/Pb, where pb is the background density. We denote the 
corresponding fields for the galaxies (or halos) by S g (x). The 



two fields are smoothed with the same window function of 
scale R s . We assume that biasing is a local process, which 
means that the galaxy density field 8 g is related to 8 within 
the local smoothing volume. 



2.1 Descriptions of Non-linear, Stochastic Biasing 

The simplest possible model for biasing is strictly linear and 
deterministic: 



S g (x) — bS(x). 



(1) 



A less restrictive definition of linear biasing follows from the 
theor y of biasing of density peaks in a Gaussian random 
field (Kaiser 1984; Bardeen et al. 1986), which predicts that 
the galaxy-galaxy and matter-matter correlation functions 
are related in the linear regime by a constant multiplicative 
factor, 



CggM = b CmmM. 



(2) 



In both cases, b is referred to as the "linear biasing param- 
eter". Eqn. g follows from Eqn. |l] but the converse is not 
true. 

Linear deterministic bias can serve only as a crude null 
hypothesis: it is not a self-consistent model, it has no the- 
oretic al motivation and it se ems inconsistent with observa- 
tions. Dekel & Lahav (1998, hereafter DL) have proposed 



a general formalism for non-linear, stochastic local biasing, 
which we shall adopt in this paper. We shall denote the 
one-point probability distribution functions (PDF) for the 
matter and galaxy density fields as P(8) and P{5 g ). By def- 
inition, both distributions have zero mean, and the corre- 
sponding variances are a 2 = (S 2 ) and a g ss (S g ). The local 
biasing relation between galaxies and matter is treated as a 
non-deterministic process, specified by the conditional bias- 
ing distribution P(S S \S). The mean biasing function b(8) is 
then defined by 



b(5)S = {8 g \S) = / dS g P(5 g \S)S, 



(3) 



This function fully characterizes the mean non-linear biasing 
and reduces naturally to the linear biasing relation Eqn. jj] 
if b is independent of 8. 

Useful statistics for characterizing the mean biasing and 
its non-linearity are the moments 



b= (b{S)8 2 )/a 2 
and 

b 2 ss (b 2 {5)5 2 }/a 2 



(4) 



(5) 



To describe the statistical nature of the biasing relation, we 
define the random biasing field 

e = fl -<*«!*>, (6) 

and its variance, the biasing scatter function 

a 2 b (S) EE {e 2 \S)/a 2 . (7) 

Averaging over 5, one obtains the biasing scatter at,: 

al ee {e 2 )/a 2 . (8) 

To second order, the three parameters b, b and <7b char- 
acterize any local, non-linear and stochastic biasing relation. 
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The parameter b is simply the slope of the linear regression 
of 5 g on 8, and as such it is a natural generalization of the 
linear biasing parameter of Eqn. [l]. The ratio b/b quanti- 
fies the non-linearity of the mean biasing relation and <Jh/b 
independently measures the scatter. Thus the non-linearity 
and stochasticity can be studied separately. The formalism 
reduces in a natural way to linear biasing when b — b and 
to deterministic biasing when <Tb = 0. We do not investigate 
moments higher than second order in this paper, but the 
formalism is readily generalizable. 

Various statistics describing biasing are used in the liter- 
ature and commonly referred to without distinction as "the" 
bias parameter. In general, these statistics are not equiva- 
lent. Here we follow the notation of DL which naturally gen- 
eralizes the linear deterministic biasing parameter of Eqn. [j] 
into the mean biasing parameter b of Eqn. ^. The biasing 
parameter relating the variances at a given smoothing scale 
is denoted b va , r , and it can be expressed in terms of the three 
basic parameters above: 



ZJL 
a 



= b 



_ + 1± 

b 2 b 2 



1/2 



(9) 



A complementary parameter is the linear correlation coeffi- 
cient 

1 + ^ . (10) 

b 2 b 2 



(S s 8) 



Note that both 6 V ar and r mix nonlinear and stochastic ef- 
fects. We denote the r atio of co rrelation functions at a given 
separation by 6 5 = i/£ gg /£ mm . 

2.2 Models of Biasing 

2.2.1 Galaxy Conserving Models 

The bias and linear correlation coefficient are expected to 
evolve with time in accord with the evolution of cluster- 
ing. If the galaxies behave as test particles in the mat- 
ter density field and their numbers and intrinsic proper- 
ties are conserved, then they satisfy the continuity equation 
The par ameters & var and r then evolve as 



( [Tegmark fc Peebles 199$ ) 

bo = [(1-D) 2 + 2D(1 
= [(1 - D) + Dbr]/b , 



D)br + D 2 b 2 } 1/2 



(11) 
(12) 



where bo and ro are the present day quantities, b and r are 
the corresponding values at some earlier redshift z, and D 
is the linear growth factor at z (normalized to unity today). 
If there is no correlation between galaxies and mass (r = 1), 
this expre ssion reduces to the model originally proposed by 
Pry (199q ). We shall refer to this type of model as "galaxy 
conserving" (GC). 



2.2.2 Hierarchical Merging Models 

Mo & White (1996, hereafter MW) developed a model for 
the biasing of virialized dark-matter halos with respect to 
the underlying matter distribution. Their model is based on 
an extension of the Press-Schechter approximation (Press & 
Schechter 1974), and accounts for halo merging and contin- 
uous formation of new halos. They derived an expression for 



Table 1. Simulation Parameters. Prom left to right: the cos- 
mological mass-density parameter Q, the Hubble constant h = 
-Ho/(100km s _1 Mpc -1 ), the linear rms mass density in a sphere 
of radius 8h~ 1 Mpc erg, the box size L, and the mass per particle, 



Model 


n 


h 


CT8 


L (h- 1 Mpc) 


m p (h- 1 M Q ) 


tCDM 


1.0 


0.5 


0.6 


85 


1.0 x 10 10 


ACDM 


0.3 


0.7 


0.90 


141 


1.4 x 10 10 



the mean nonlinear biasing function b(5; M, z, R a ) of halos 
with mass M at a redshift z for a smoothing scale R s (MW 
eqn. 19-20). In the vicinity of | <5| ^> 5i, where Si = S c (l + z) 
is the density of a collapsed object, this reduces to the simple 
linear biasing factor: 



b(M,z) 



1 + 



Si 



(13) 



where v = 8i/a(M), and o(M) is the rms mass variance. 
MW showed that their model predictions for these quanti- 
ties are in good agreement with the results of N-body sim- 
ulations with scale-free initial power spectra, P(k) oc k n . 



2.2.3 Shot Noise and Stochasticity 

In empirical measures of the scatter in the biasing relation, 
shot noise is inevitably mixed with the physical sources of 
stochasticity. Removing this shot noise is not straightfor- 
ward because halos or galaxies correspond to rare peaks in 
the density distribution, and their selection is far from a 
Poisson process. Let us imagine however a deterministic bi- 
asing scheme with a known form (<5g|<5) (for example, mea- 
sured from a simulation) . The expected variance in the num- 
ber of galaxies A in a smoothing window with volume V 
and overdensity S is ^ 2 = {(N - (A)) 2 ) = A = nV(l + 5 S ), 
where n is the global average number density of galaxies. 
The variance of the galaxy field 5 S as a function of S is then 
0-9 (*) = /(nV) 2 = (l + $g(5))/(nV). We shall refer to this 
simple model as the "conditional Poisson" model. 



3 SIMULATIONS 

The dissipationless iV-body simulations and the semi- 
analytic method used to model galaxies within these sim - 
ulations are described in detail in Kauffmann et al. (1998a) 
Here, we briefly summarize only their main features. 



3.1 A^-body Simulations 

Special A-body simulations were run for this project 
(termed "GIF") using the version of the adaptive particle- 
particle particle-mesh (AP 3 M) Hydra code developed as 
part of the VIRGO supercomputing project. The simula- 
tions have A p = 256 3 particles and 512 cells on a side, 
and a gravitational softening length of 30/i _1 kpc (h = 
^/(lOOkms" 1 Mpc^ 1 )). Here we analyze only two of the 
four cosmological models simulated for the GIF project (see 
Table §). One model, rCDM, has Q = 1 and the other, 
ACDM, has SI = 0.3 and a cosmological constant Q,a = 0.7. 
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The initial power spectra of fluctuations both have a shape 
parameter F = 0.21 and are normalized to approximately 
reprod uce the number density of clusters at the present 



The star formation history of each galaxy is convolved 
with stellar population synthesis models to obtain total lu 



minosities in any desired filter band. The models of Bruzual 



epoch. Dark matter halos were identified using a standard & Chariot (1998) are used, assuming that all stars have so- 



friends-of- friends algorithm with a linking length 0.2 times 
the mean interparticle separation, thus corresponding to a 
density contrast of ~ 125 at the outer parts of the halos. 



We do :iot include any halos smaller than 10 particles in our 



analysis, as tests show that these halos are not stable over 
many output times. 



3.2 Semi- Analytic Modelling of Galaxies 

A "merger tree" is constructed for each halo identified at z = 
by searching for its progenitor halos at earlier redshifts. A 
halo at an early redshift Zi is defined to be a progenitor of 
a halo at a later redshift zo < Z\ if more than half of the 
particles of the progenitor and its most-bound particle are 
included in the halo at Zq. The progenitors are the halos that 
will merge together to form the parent halo at z — 0. If a halo 
is at the top level of the hierarchy (i.e. it has no progenitor 
with at least 10 particles), then it is assumed to contain hot 
gas at t he virial temperature of the halo, This gas is allowed 



to COol larrnrrling tr. tViP r-Hrliati-irp pnnling timpgpalp anrl thp 

cold gas is assumed to settle into a galactic disk and form 
stars. 

Th e atar formation rate ia given by the cxprcaaion 



77gog 

^dyn 



(14) 



where a is a free parameter, m co id is the mass of cold gas 
in the disk, and td yn is the dynamical time of the disk. The 
dynamical time is approximated as £ dy n = 0.1r v i r /V" c , where 
r v i r is the virial radius and V c is the circular velocity of 
the halo at r v i r . Cold gas may be reheated by supernovae 
feedback, with the rate of reheating given by the expression 



(15) 



where e is a free parameter, and Vo is a scaling constant. 
Reheated gas is removed from the cold gas reservoir. It may 
then be retained in the halo where it will cool again on 
a relatively short time scale, or ejected from the halo. If 
ejected, the gas is returned to the halo after the mass of the 
halo has grown by a factor of two or more. 

The N-body simulations used here do not have sufficient 
resolution to resolve sub-halos once they are incorporated 
into larger halos. Therefore, the merging rate of galaxies 
within the halos is modeled semi-analytically. When halos 
merge, the central galaxy of the largest progenitor halo be- 
comes the central galaxy of the new halo and all other galax- 
ies become "satellite" galaxies. A satellite galaxy is assumed 
to merge with the central g alaxy once its dynamical fr iction 
timescale has elapsed (see Binney & Tremaine 1987). The 



morphologies of the galaxies are determined by assuming 
that major mergers (m 8a t/TO C entral > 0.3) result in destruc- 
tion of the disk and consumption of all remaining gas in a 
starburst. The remnant is assumed to be a bare spheroid. 
Subsequent gas cooling may result in the formation of a new 
disk. The bulge-to-disk ratio at each output redshift is then 
used to divide the galaxies into rough morphological cate- 
gories. 



l ar metallici ty, and a standard Scalo initial mass function 
( [gcalo 1986 ). The effects of dust extinction on the galaxy 
spectra were investigated using the empirical recipe of Wang 



Heckman (1996) (WH), in which the face-on optical depth 



in the B-band is parameterized in the form: 

TB OC T B (L B /L* B ) , 



(16) 



with P = 0.5, r B = 0.8, and L B = 1.3x W 10 h' 2 L Q . This can 
be extended to other wavebands using a standard Galactic 
extinction curve. The extinction is then computed by assign- 
ing a random incli nation to each galaxy and using a standard 
"slab" model (see Kauffmann et aL (1998a) for details). 

The two free parameters in the galaxy formation mod- 
els (a in Eqn. [l4| and e in Eqn. [l5|) are set by requiring the 
central galaxy in a halo with V c = 220 km s _1 at the present 
epoch to have an I-band magnitude of Mr — 5 log h ~ —22.1 
and a cold gas mass of m co id ~ 10 10 Mq. This ensures that 
the zero-point of the Tully-Fisher relation for spiral galax- 



ies is in agreement with observations (e.g. Giovanelli et al 



199J ). 

Various properties of the galaxies produced using these 
techniques are summariz ed in the companion pape rs based 
on the GIF simulations, 



Kauffmann et al. (1998a ) 



maim et al. (1998b ), and Diaferio et al. (1998 ). We note 



Kauff- 



that throughout this paper, we include in our analysis only 
galaxies that reside in halos of least 10 particles. 



4 FROM HALOS TO GALAXIES 

The study of biasing is in some sense the attempt to under- 
stand the connection between the total mass density, which 
we believe to be dominated by dark matter, and the density 
of luminous galaxies. A fairly robust component of modern 
galaxy formation theory is that galaxies form within col- 
lapsed, virialized dark-matter halos. The clustering proper- 
ties of these halos relative to the underlying mass distribu- 
tion are straightforward to compute from cosmological N- 
body simulations such as the ones that we have described, 
subject only to the numerical limitations of the simulation 
at hand. Predicting how these halos are connected to visible 
galaxies is far more difficult. We have adopted the particular 
recipes described above, but other equally plausible recipes 
could lead to different results. In this section we quantify 
the resulting statistical relation between halos and galaxies 
in our simulations, and in the next section we present in 
parallel the results for the biasing properties of halos and 
galaxies. In this way we hope to gain some insight into the 
extent to which the biasing of galaxies is determined by the 
gravitational process of halo formation and to what extent 
by other astrophysical processes. 

In general, a single dark-matter halo in our simulations 
may contain several galaxies with a range of luminosities. 
The largest collapsed halos in our simulations have virial 
radii of 1.5 — 2h~ J Mpc. Therefore, on the smoothing scales 
that investigate here (R s > 4/i -1 Mpc), the clustering of 
galaxies is mainly determined by the clustering of the ha- 
los in which they dwell. Figure hi shows the joint probabil- 
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ity function of absolute magnitude and host halo mass for 
the galaxies in our simulations. The pronounced diagonal 
"ridge" is populated by galaxies that are the central object 
in their halo, and illustrates the fairly tight relationship be- 
tween luminosity and halo mass or circular velocity (i.e. the 
Tully-Fisher relation) obeyed by these galaxies. This ridge 
corresponds to a constant halo mass-to- light ratio, where the 
value of the constant is constrained to agree with the zero- 
point of the observed Tully-Fisher relation by construction 
(see Section 3.2). The galaxies that lie far off of the ridge are 
mainly satellites. Thus we see that galaxies with a given lu- 
minosity live in a wide range of environments. The differing 
appearance of the contours in the two different cosmologies 
(for example the much broader ridge in the rCDM simula- 
tions) are a direct consequence of the details of the star for- 
mation and supernovae feedback recipes, which were chosen 
to give the best agreement with observations (Kauffmann 
et al. 1 )98a). Put another way, these two cosmologies have 
different halo mass functions, so that in order for both to 
reproduce the observed galaxy luminosity function, the re- 
lationship between halos and galaxiesmust be different. The 
right-hand panels show the results of including dust extinc- 
tion. We see that this effect can substantially modify the 
relationship between mass and light in the B-band. We re- 
turn to this point later. 

Figure ti shows the mass-to-light ratio of halos in our 
simulations in the B and I bands. The effective mass-to- 
light ratio when dust extinction is included is also shown. 
The mass-to-light ratio varies by more than one order of 
magnitude from the smallest to largest halos in the simu- 
lations, and has a large scatter at a fixed halo mass. The 
characteristic mass M», defined as the mass that is just 
becoming non-linear at a particular epoch (more precisely 
ff(M»(z)) = 5 c (z)), is M» = 2.9x10 
model and Af, = 1.5 x W^h' 1 M { 

(see figure []). The exponential turn-over in the halo mass 
function thus occurs at a mass much larger than that of halos 
that typically host L* galaxies. Therefore the sharp increase 
in M/L at ~ 10 13 /i -1 Mq is necessary in order to produce 
a luminosity function with a "knee" at L„. This is generally 
believed to be the result of inefficient cooling in large ha- 
los, and in these simulations it is achieved by turning off gas 



+2 



h,- 1 Mq fortherCDM 
for the ACDM model 



Kauffmann et al 



cooling in halos larger than 350 km s _1 (see 
1998ahJ We may also note that the difference between the 
mean B-band and Fband mass-to-light ratio increases with 
halo mass, indicating that larger halos tend to host galaxies 
with redder colours. 

The distribution P(n\Mh) of the number of bright 
galaxies per halo as a function of halo mass is shown 
in figure W. Halos with masses between 10 12 /i -1 Mq and 
10 13 /i _1 Mq contain one bright galaxy on average, with 
a relatively small dispersion. Larger mass halos (> 5 x 
10 13 fo -1 Mq) typically contain more than one galaxy per 
halo, and the number of galaxies per halo has a larger dis- 
persion; these halos represent groups or clusters of galaxies. 



5 BIASING OF HALOS AND GALAXIES 

In this section we present our main results concerning the 
biasing relation between the overdensity fields of halos, or 
galaxies, and that of the underlying dark matter. 



Figure |H| shows the joint distribution P(8 g ,8) for ha- 
los, and figure ^ shows the same for galaxies. All the den- 
sity fields have been smoothed with a top-hat filter of ra- 
dius 8 h~ x Mpc (hereafter T8). We have selected halos with 
masses larger than 1.0 x 10 12 /i _1 Mq, and galaxies with 
Mb — 5 log h < —19.5, at several redshifts as indicated on the 
figures. The conditional mean and standard deviation func- 
tions are shown, corresponding to b(8) and Ob(S) of Eqn. ^ 
and [5]. Shown for reference is also the linear biasing approx- 
imation <5 g = bS, with b — b. A log-log density plot is used in 
order to stress the behavior in the regions of under-density. 
Note that in this specific presentation, linear biasing appears 
as a curved line. 

The results for halos and galaxies are qualitatively sim- 
ilar. This is somewhat surprising given the rather complex 
relationship between the two that we saw in the previous 
section. However, recall from figure ^ that on average there 
is approximately one galaxy brighter than —19.5 + 5\ogh 
per ~ 10 12 fo _1 Mq halo. These halos are much more nu- 
merous than the larger mass halos that host larger numbers 
of galaxies, so they tend to dominate the appearance of a 
number-weighted joint probability such as the quantity rep- 
resented in these figures. We also note that the contours for 
the two cosmological models differ in their details but overall 
they appear qualitatively similar. 

The favored environment for halo/galaxy formation, 
namely, the peak in the halo/galaxy overdensity distribu- 
tion, occurs in regions where the matter density is close to 
its mean value (1 + 5 = 1). In the vicinity of 1 + 6 ~ 1, the 
overdensity of halos/galaxies follows that of the underlying 
matter, b(8) ~ 1. But, in general, the linear biasing approx- 
imation is not an accurate description of the mean biasing 
relation. 

The mean biasing function shows a robust characteris- 
tic behavior in the under-dense regions (1 + 8 < 1); the local 
slope of b(S) is steeper than the global slope b, and is larger 
than unity even when b < 1. This leads to voids empty of 
galaxies below some finite matter underdensity (1 + 8 > 0); 
namely, the galaxy-formation efficiency drops to zero be- 
low a certain mass-density threshold. We can say that the 
galaxies in the voids are always "positively biased", in the 
sense that locally b > 1, while they can be either biased or 
"anti-biased" (6 < 1) in the high-density regions. 

In the over-dense regions (1 + 5 > 1), the behavior 
changes with redshift. At z = 0, the local slope of b(S) is 
smaller than unity, driving b to below unity. At higher red- 
shifts the local slope becomes larger, driving b to values as 
large as 3 to 7 by z — 3. The redshift dependence of the 
biasing relation is striking; we return to it in 5.4. 

The stochasticity in the biasing scheme is evident from 
the spread in S g at fixed 8. We will discuss it further in the 
next section. 



5.1 Mass and Luminosity Dependence 

In this section we investigate the dependence of biasing on 
halo mass and galaxy luminosity, for a fixed smoothing win- 
dow (T8) and at the present epoch (z = 0). Figure |^ (bot- 
tom) shows the mean biasing function for halos selected with 
different mass thresholds, in the rCDM simulation. We see 
that massive halos are more biased; the curves all cross at 
8 ~ 0, massive halos have higher overdensities than smaller 
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mass halos in regions with S > 0, and smaller mass halos 
have higher overdensities where 8 < 0. 

The top panel shows the variance of the conditional 
biasing distribution {e 2 1 <5) minus the mean shot noise cor- 
rection l/(nV s ), where n = N/Vbok. is the average number 
density of halos/galaxies above the given threshold in the 
simulation, and V s is the volume of the smoothing window . 
The conditional Poisson model discussed in Section 2.2.3 is 
shown for comparison. In regions of lower than average over- 
density, the scatter is generally smaller than the mean shot 
noise, and in overdense regions it is larger. 

This can be understood by re ferring to the expression 
for the variance of counts in cells (Peebles 1980): 



^ 2 = ((N-nV) 2 } =nV + n 2 / dV 1 dV 2 .£, 



(17) 



The variance of the conditional biasing distribution is then 
Hz/inV) 2 . This reduces to the usual shot noise value 1/V nV 
in the absence of correlations (£ = 0), but it leads to a 
scatter larger than the mean shot noise when the integral 
over £ is positive, and smaller than the mean shot noise 
when it is negative. In underdense regions (S < 0), halos are 
anti-correlated (£ < 0) and the volume- averaged correlation 
function, £, is negative. In overdense regions (5 > 0), halos 
are positively correlated and £ is positive. 

The corresponding results for the ACDM model are ex- 
tremely similar (and are therefore not shown) because the 
mass dependence enters only through the shape of the power 
spectrum, which is similar by design for the rCDM and 
ACDM simulations. 

Figure ^ summarizes the mean and scatter of the condi- 
tional biasing relation for galaxies with different luminosity 
thresholds. Two notable features are evident. First, we see 
no significant change in the biasing relation for galaxies with 
absolut e- magnitude thresholds in the range —18.4 + 51ogfe 



halos and is nearly constant over the mass range. The non- 
linearity is even smaller for galaxies, and shows a trend with 
galaxy luminosity. The stochasticity, <Jb/b, shows a modest 
trend, increasing for higher mass halos or brighter galax- 
ies. As remarked before, the stochasticity is actually larger 
for halos than for galaxies. The linear correlation coefficient 
r = 6/&var, which, like & va r, mixes nonlinear and stochas- 
tic effects, is roughly constant over these ranges of mass or 
luminosity and has a value r ~ 0.9. 

We further examine the clustering of halos of different 
masses compared to the dark matter in figure jlo| which 
shows the auto-correlation functions of halos and dark mat- 
ter in the rCDM simulations. More massive halos show cor- 
relations of higher amplitude and therefore more positive 
biasing. The turn-over of the halo correlation functions at 
small radii is caused by exclusion effects due to the finite spa- 
tial extent of the halos. When using conventional halo finders 
based on friends-of-friends or spherical over-density, halos 
that lie too close together are "merged" into a single halo. 
The spatial scale at which this effect becomes important 
depends on the mass of the halo, according to the relation 
between virial radius and mass (M v ir oc r^ iT ). Thus on small 
scales, one expects anti-correlations due to this exclusion ef- 
fect to lead to negative £ in regions of high halo density, and 
thus to a reduction in the scatter, but this effect is negligible 
for the smoothing scales considered here (R s > 4/i _1 Mpc)^. 
Note that the fact that the observed galaxy correlation func- 
tion is a nearly perfect unbroken power-law on sub-Mpc 
scales indicates that there must be a significant contribu- 
tion from galaxy pairs that dwell within a common halo. 

The correlation function of galaxies with different lu- 
minosity thresholds is sh own in figure As noted by 



Kauffmann et al. (1998a), the clustering amplitude of the 



galaxies in the simulations is nearly independent of the lu- 



minosity threshold, in contrast to the findings of Willmer 



(well b elow L t ) to —19.9 + 5 log h (above L*). We return to 
this point in a moment, when we discuss figure [nt Second, 
the scatter compared to the average shot noise is comparable 
to the corresponding quantity for halos. On the face of it, one 
might think that the complex physics of gas cooling, star for- 
mation, merging, supernovae feedback, etc, even as simply 
modelled in our simulations, would lead to a larger scatter. 
However, on the other hand, the nature of the physics asso- 
ciated with these processes might actually lead to a stronger 
correlation of luminosity with the local matter density than 
that shown by the halos. For example, gas can only cool and 
form stars in regions with sufficiently high density, the star 
formation efficiency is assumed to be proportional to the 
density, supernovae feedback is less efficient in halos with 
deep potential wells (i.e. high density), and starbursts oc- 
cur preferentially in high density regions. While these lines 
of argument appear plausible, more detailed investigation of 
the physical source of scatter in the biasing relation that we 
obtain both for galaxies and halos is clearly in order. 

The three characteristic parameters of the biasing 
scheme (see Section 2.1), measuring mean biasing, nonlin- 
earity, and stochasticity, are shown in figure ^ as a function 
of halo-mass threshold or galaxy luminosity threshold. The 
mean biasing is characterized by 6, and also shown are b val 
and 65. The fact that & var tends to be systematically larger 
is due to the presence of nonlinearity and stochasticity. The 
nonlinearity, as measured by b/b, is less than 10 percent for 



Da Costa, fc Pellegrini (1998 ), who observed pronounced 



luminosity-dependent bias in the SSRS2 redshift survey 
They found that brighter galaxies are more strongly clus- 
tered than fainter galaxies, and thus have a longer correla- 
tion length ro as shown on the figure. 

Our result is not surprising given the weak mass de- 
pendence of the biasing of galaxy-mass halos (figure [J). 
However, this introduces an apparent problem when com- 
pared to observations. The difference between the faintest 
and brightest magni tude thresholds for the observed galaxy 
samples analyzed by Willmer, Da Costa, & Pellegrini (1995) 
is about 1.5 magnitudes. This corresponds to a factor of four 
in mass if the mass-to-light ratio is constant and its value is 
dictated by the observed Tully-Fisher relation. On the other 
hand, we saw from figure ^ that there is only a ~ 10% in- 
crease in the mean biasing as a function of halo mass over 
this range, and even this weak mass dependence is diluted 



* The exclusion effect results in a reduction of the scatter by a 
factor that scales like (r ex / iJ s ) 3 , where r ex is the scale on which 
the exclusion effects are important. Here (r ex /ij s ) 3 ~ 2 X 10 — 2 — 
2 X 10~ 3 depending on the halo mass, for R B = 4/1" 1 Mpc. 
t The galaxy catalogs used here are not identical to the ones used 
in Kauffmann et al. (1998g] ), which explains the small differences 
in the correlation functio ns shown here and Figure 10 and 11 of 
Kauffmann et al. (1998a). 
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by the presence of low-luminosity satellite galaxies in large- 
mass halos (see figure Q), yielding no luminosity dependence 
in our simulated galaxy biasing. 

A clue for a resolution of this problem may come from 
the fact that the mass dependence of halo biasing becomes 
much stronger in the regime where M > M» (figure []) . If we 
increase the mass-to-light ratio, pushing the bright galaxies 
into larger mass halos, we would obtai n stronger luminosity 



depend ence, as wa s in fact found by Kauffmann, Nusser 



& Steirmetz (1997) when they effectively normalized their 
models with a higher mass-to-light ratio. Unfortunately, this 
would lead to an apparent contradiction with the observed 
Tully-Fisher relation. 

A more promising way to achieve the desired effect is 
by inclusion of dust in our considerations. We note that the 
observed Tully-Fisher relation has been corrected for dust 
extinction using observationally determined inclinations for 
each galaxy. These corrections can be quite large in the B- 
band, where most of the samples used to study galaxy clus- 
tering are selected. Moreover, there is evid ence that more lu 



minous galaxies suffer la rger extinctions (Fully et al. 1999 
Heckman 199^ ). Studies of galaxy clustering do 



Wang 

not include corrections for dust extinction because inclina- 



tion estimates are generally not available. Kauffmann et al 
(1998a| showed that including dust extinction can signifi- 
cantly modify the galaxy-galaxy correlation function found 
in these simulations. 

The tCDM model suffered from an excess of bright 
galaxies compared to the observed B-ban d luminosity func- 



tion when dust extinction was neglected. Kauffmann et al 
(1998a| found that including dust extin ction using the em- 



pirical recipe of |Wang fc Heckman (199tj ) led to an improved 
luminosity function, but still with an excess, especially on 
the bright end. We now tune the parameters of the WH 
recipe in order to obtain the best possible fit to the ob- 
served luminosity function. For the rCDM model, we find 
that using (3 = 0.2 and setting t b = 2.0 in eqn. |l^ gives the 
best results. For the ACDM model, the luminosity function 
already has a small deficit of galaxies compared to obser- 
vations even without any dust correction. Despite this, we 
show the results of applying a dust correction to the ACDM 
models using the fiducial values of the WH parameters (as 



in Kauffmann et al. 1998a) 



The inclusion of dust extinction changes the effective 
mass-to-light ratio of the halos in the simulations (see fig- 
ure |l|), in the desired sense: galaxies of a given observed 
luminosity now dwell in larger mass halos. We see in the 
right-hand panels of figure hll that this leads to a luminosity 
dependence that is qualitatively similar to the observed de- 
pendence, especially in the rCDM case (unfortunately, the 
number of bright galaxies in our simulation box also becomes 
quite small, and the correlation function becomes rather 
noisy, but the trend is clear). Note also that although the 
inclusion of dust increases the correlation amplitude of the 
galaxies in the rCDM model, bringing the results into better 
agreement with the observed galaxy correlation function, the 
clustering amplitude for bright galaxies is still not as high 
as the SSRS2 observations. Conversely, with no dust cor- 
rection, the clustering amplitude of galaxies for the ACDM 
model is already comparable to the observational results, 
and inclusion of dust hardly makes a difference. It may be 
that a cosmology with an intermediate value of fio — 0.5 will 



give the best overall results when realistic extinction due to 
dust is included. However, a much larger correction seems 
to be required in order to match the correlation function 
of very bright galaxies as provided by these observations. 
These results highlight the importance of obtaining large 
redshift surveys selected in near IR bands in order to reduce 
the sensitivity to dust extinction. 



5.2 Type Dependence 

We study the type dependence of biasing with T8 smooth- 
ing, for galaxies brighter than Mb — 51og/i = —18.4, at 
z — 0. Here, instead of quantifying the bias of galaxies with 
respect to mass as in the previous section, we show the bias- 
ing of different types of galaxies relative to the overall galaxy 
population. This relative bias is of particular practical inter- 
est because it can be compared directly with the results of 
observations. 

We can classify different types of galaxies in our simula- 
tions in several ways. For example, we identify "early" and 
"late" types according to the bulge-to-total luminosity ratio; 
galaxies of Lbui gc /£tot > 0.4 are identified with early type 
galaxies (Hu bble type E-SO), and th e rest with late types 



Kauffmann et al. (1998a). Similarly, we can di- 



(S-Irr), 

vide the galaxies in terms of colour. Here, we classify galaxies 
with B — V > 0.8 as "red" and the remainder as "blue". Fig- 
ure ^| (top) shows the joint density distribution for galaxies 
of early and late types relative to the galaxy population as a 
whole at the same magnitude limit. Early type galaxies are 
biased (b > 1) compared to the overall population, wheras 
late type galaxies are slightly anti-biased (b < 1) We ob- 
tain comparable results, with stronger bias and anti-bias, 
when galaxies are divided in terms of their colors, as shown 
in the bottom panels. The results shown are for the rCDM 
simulation, and are similar for the ACDM simulation. 

Once again, we can understand this result in terms of 
the masses of the dark matter halos in which these different 
types of galaxies are found. Figure ^| shows the distribu- 
tions of host halo masses for early and late type galaxies 
and red and blue galaxies. Although each type occupies a 
broad range of halo masses, the mean halo mass occupied 
by early/red types is significantly larger than that occu- 
pied by late/blue types. Several physical effects included 
in the semi-analytic modelling may contribute to this re- 
sult. Larger mass halos are more likely to have suffered a 
recent major merger, which is assumed to destroy the disk 
and create a bulge. In addition, gas cooling ceases in large 
mass halos, cutting off the supply of new gas and subse- 
quently the star formation. Therefore the galaxies in these 
halos redden, are not able to form new disks and so remain 
bulge-dominated. Note that the separation between the dis- 
tributions of red and blue galaxies is larger than the separa- 
tion between early and late types (figure ^) , explaining why 
the relative biasing is stronger for the former (figure p^ ). 

The mean relative biasing parameters found in the sim- 
ulations are summarized in Table ^. These parameters are 

where 



2.1 



the equivalent of b and 6 var as defined in Section 
the matter field 5 is replaced with the field of late/blue type 
galaxies and <5 g is replaced with the field of early/red type 
galaxies. The difference between & rel and b^ar again indicates 
the presence of nonlinearity and stochasticity in the relative 
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Table 2. Relative bias of early to late and red to blue galaxies 
in the simulations on a scale of 8 h~ 1 kpc. For each parameter, 
the first value is for tCDM and the second is for ACDM. 



type 


6 rcl 


6 rel 

var 


r 


-el 


early/late 


1.4 1.1 


1.5 1.3 


0.93 


0.87 


red/blue 


1.6 1.6 


1.8 1.8 


0.87 


0.87 



biasing relation, and is also reflected in the relative linear 
covariance parameter r rcl . 

Observationally, the stronger clustering of early type 



galaxies is well known (Dressier 1980 ; Hermit et al. 1996 



Willmer^ Da Costa, fc Pellegrini 1998; Tegmark fc Brom- 



ley 199|j|). In order to compare with theoretical predictions, 
care must be taken to account properly for the transfor- 
mation from redshift space to real space, as early and late 
type galaxies are known to be affected differently by red- 
shift distortions. In real space, Willmer, Da Costa, fc Pcl- 
legrini 1998) find & C ariy/&iatc = 1. 18 ± 0.15, Lovcday et al 
(1995^ ind 6 car i y /6i atc = 1.33 and GJuzzo et al. (1997) find 
focariy/^iato = 1.68. Not surprisingly, the range of observa- 
tional values seems to depend on the importance of clusters 
in the sample, with higher values obtained in samples that 
include many rich clusters. The quoted observational statis- 
tics correspond to 6 V ar- The values we obtain (1.3-1.5) are 
in reasonable agreement with the observational values. 

It is also fairly well established that red galaxies are 



general ly more clustered than the blue population (Landy. 
Szalay, & Koo 1996; Tucker et al. 1997), but a quantitative 



comparison is more difficult, because most of the values in 
the literature are obtained from the angular or redshift space 
correlation functions, and a variety of c olour bands are used. 
Willmer, Da Costa, & Pellegrini (1996) quote a relative vari- 
ance bias in real space bred/6blue 1.40 ±0.33 at 8h~ Mpc, 
and suggest that there is evidence for scale dependence in 
the relative bias. However, they used a colour threshold of 
B — R = 1.3 so it is not directly comparable with our results. 

In principle, the type dependence could be further stud- 
ied as a function of redshift, scale, luminosity, or environ- 
ment. The samples currently available from both simula- 
tions and observations are too small to obtain proper statis- 
tics after this sort of subdivision, but this may be a pow- 
erful discriminatory tool for galaxy formation models once 
larger simulations and larger redshift surveys are available. 
Another promising approach would be to categorize galax- 
ies using automatic spectral cl assification methods such as 
Princip al Component An alysis (Connolly et al. 1995; Folkes, 
Lahav, & Maddox 1996) and study the corresponding rela- 
ti ve bias. An example of a co mparison of this sort appears 
in [Tegmark & Bromley (19981). 



5.3 Scale Dependence 

Figure [w] and |l5| show the scale dependence of biasing for 
a fixed halo mass threshold (M > 10 12 /i _1 A/0) or galaxy 
luminosity (Mb — 51og/i < —19.5) at the present epoch 
(z = 0). The results shown are for top-hat smoothings of 4, 
8 and 12h~ Mpc. This figure shows the mean and scatter 
of the conditional biasing relation for halos in the rCDM 
simulation (again, the results of the ACDM simulation are 



nearly identical and are not shown). The scale dependence 
of the mean biasing is weak over this range, while the scat- 
ter shows an expected scale dependence. Once again, the 
results for galaxies are similar to the results for halos. This 
suggests that, at least in our simulations, the scale depen- 
dence (or lack thereof) of biasing of bright galaxies is mainly 
determined by the gravitational physics of halo formation. 

Figure [l6] presents the scale dependence of the bias- 
ing parameters. Here we see that the mean biasing actually 
increases slightly with smoothing scale for both halos and 
galaxies. As expected, the non-linearity and stochasticity 
both decrease with increasing smoothing scale, such that 
the biasing relation converges towards linear deterministic 
biasing for large smoothing scales. 

5.4 Redshift Dependence 

We now return to the redshift evolution of biasing, which 
is the most striking feature in figure JE] and and is of 
particular interest given the recent observational develop- 
ments at high redshift. We use a fixed halo-mass threshold 
(M > 10 12 /i _1 M Q ) or luminosity threshold (M s -51ogft < 
— 19.5) and smoothing scale (T8). We already noted that the 
shape of the biasing relation changes noticably, and in par- 
ticular the slope b increases dramatically with redshift. The 
rapid evolution towards higher bias at high redshift has been 
noted in several recent works (B agla 1 9981; Wechsler et al 



1996 



Colin et al. 1998 



Katz, Hernquist, & Weinberg 1996 
Bfanton^cT*aL*F999[ ). Here we discuss several other interest- 
ing properties of the redshift evolution of biasing of both 
halos and galaxies. 

The tCDM and ACDM models, which showed very sim- 
ilar mass and scale dependences at z = 0, show quite differ- 
ent redshift evolution. The evolution is more pronounced in 
the Q = 1 rCDM model than in the low-fi ACDM model, 
as expected. This is because the mass and scale dependence 
at fixed redshift are determined mainly by the shape of the 
power spectrum, which is deliberately similar for these two 
models, while the redshift evolution depends on the growth 
rate of clustering, which is quite different in these two mod- 
els. 

In addition, we can see in figure ^ and ^ that the evolu- 
tion of the biasing relation for galaxies is weaker than that 
for halos. This is because of the redshift dependence of the 
efficiency of star formation and hence of the relationship be- 
tween halo mass and galaxy luminosity. Halos are typically 
denser at high redshift, and in these models we have effec- 
tively assumed that the star formation rate is proportional 
to the averag e halo density (see figure 8 and 9 of Kauffmann 
et al. (1998b| ) and Eqn. hj of this paper). Therefore, when 



we select galaxies with a fixed luminosity threshold, at high 
redshift we include galaxies residing in smaller mass halos 
(and thus of lower biasing). Note that this effect also re- 
duces the difference in the evolution of biasing within the 
two cosmological models, i.e., the additional physics asso- 
ciated with star formation conspires to make the redshift 
evolution of biasing less discriminatory to cosmology. 

Figure [l?] summarizes the redshift evolution of the mean 
bias, nonlinearity, and stochasticity parameters. Note the 
increasing difference between the mean biasing statistics, 6, 
&var and &£, with increasing redshift, for the rCDM model; 
b differs from b% by a factor of two at z = 3. This differ- 
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ence reflects the increase in stochasticity and non-linearity 
with redshift in the rCDM model, as seen in the lower pan- 
els of the figure. The increase in stochasticity in this model 
is dominated by the increase in shot noise due to the de- 
crease in number density of halos/galaxies between z = 1 
and 3. The difference in and 6 V ar is due to the stronger 
scale dependence at z — 3 (b$ reflects the biasing at a fixed 
scale, whereas 6 var is an average over different scales). In 
the ACDM simulations, the non-linearity and stochasticity 
instead decrease with redshift, and therefore b, 6 var and foj 
agree within 30 percent at z = 3. Thus, the evolution of the 
biasing parameters is characterized in the two models by a 
minimum in the nonlinearity and stochasticity, which occurs 
at z ~ 1 in the rCDM cosmology and z ~ 3 in the ACDM 
cosmology. 

Some insight into the origin of this minimum may be 
gained by examining the time evolution of the correlation 
functions of dark matter, halos, and galaxies from z — 3 
to z = 0, shown in figure [l^. The clustering amplitude of 
the dark matter increases monotonically as time progresses. 
However, halos of a fixed physical mass correspond to rarer 
peaks in the density field at higher redshift, so the clustering 
amplitude of these objects decreases monotonically as time 
moves forward. The changing mass scale corresponding to 
galaxies of a fixed B-band luminosity discussed above leads 
to a rather different behaviour for the magnitude-limited 
galaxies. In both models, the clustering amplitude decreases 
as we move backwards in time from z = to z = 1. It 
starts increasing at a redshift of about z — 1 for rCDM and 
z — 3 for ACDM. This effect is discussed in much more de- 
tail in (Kauffmann et al. (1998h| ), who found that the magni- 
tude of this "dip" in the clustering amplitude or correlation 
length depends on the way in which galaxies are selected 
as well as on the cosmology. In both models investigated, 
the redshift at which the the correlation amplitude begins 
to rise corresponds approximately with the minimum in the 
non-linearity and stochasticity of the biasing relation noted 
above. This correspondance is intriguing but its significance 
is not immediately obvious. 



6 COMPARISON WITH ANALYTIC MODELS 
OF BIASING 

Analytic models of biasing provide an important counter- 
part to A r -body simulations, which have limited resolution 
and volume. In this section we evaluate several biasing mod- 
els by comparing them with the results of our simulations, 
focussing espe cially on the redshift dependence. Numerous 
recent papers ( Matattese et al. 1997^Moscarfin^ct^L_^998l 



Adelbergcr et al. 199$ poles et al. 1998| ; Mo, Mao, fc White 
199S ; MaglTo7clietTi^taL*r999[ ) have made use of versions 
of these analytic models to make predictions about galaxy 
clustering at high redshift and to interpret high-redshift ob- 
servations. 

We first in vestig ate the Mo & White m odel, summa- 
rized in Section 2.2.2| . In their original paper, Mo & White 



(1996) showed that their model agreed well with the results 



of simulations with scale-free power spectra. It is useful to 
investigate their model in the context of a more realistic 
power spectrum, and to show results in terms of the physi- 
cal mass scales and redshifts that we expect to correspond 



to observable galaxies. For technical convenience, we dis- 
play the biasing relation as a function of mass threshold, 
while MW originally provided predictions for halos of a spe- 
cific mass. To compute the MW prediction for all masses 
above a threshold, we simply perform a weighted average 
over the original expression (MW eqn. 19) , using the Press- 
Schec hter formula for the num ber density of halos of a given 
mass (Press & Schechter 1974). 

Figure [H^ shows the comparison of the MW predictions 
for the mean biasing relation with the simulation results 
(rCDM) for halos at different mass thresholds, smoothing 
scales, and redshifts. Apparently the MW model does well 
at reproducing the changing shape of the biasing relation 
with mass, scale, and to a lesser extent, redshift. At z — 0, 
the MW model tends to slightly underpredict the halo over- 
density in regions of very low matter overdensity. At larger 
redshifts, the MW model fails to follow the rapid growth 
of the slope b with redshift. The progressively larger de- 
viation of the MW predictions from the simulation results 
with increasing redshift is probably the result of inaccura- 
cies in the underlying Press-Schechter and extended Press- 
Schechter approximat ions, which are known t o have a similar 
redshift dependence ( [Somerville et al. 199£ ). However, the 
MW model does correctly predict the qualitative features of 
the evolution, including both the change in shape and the 
progression towards higher biasing. 



As we mentioned in Section 2.2.2, when the overdensity 



within the smoothing window is low compared to the criti- 
cal overdensity for collapse at a particular epoch, the MW 
model reduces to a very simple expression for the linear bias 
factor b(M, z) as a function of halo mass and redshift (MW 
eqn. 20, our eqn. jl3|). Figure ^o] shows the mass dependence 
of the halo biasing parameter fe V ar at several redshifts for 
the simplified MW model and the rCDM simulations. Even 
this simple expression agrees remarkably well with the mass 
and redshift dependence of biasing of halos in our simula- 
tions. The mass range of halos that can be usefully studied 
in our simulations is limited on the low mass end by our 
resolution (M rcs = 2.0 x 10 11 Mq) and on the high mass end 
by shot noise. In the figure, we show the simulation results 
both with and without a standard Poisson correction for 
shot noise. Particularly for high mass halos, the sampling 
is far from Poisson and the shot noise correction is likely 
to be inaccurate. Another interesting thing to note is that 
the mass dependence of biasing over scales comparable to 
the size of galactic halos (~ 10 11 — 10 12 /i -1 Mq) is much 
stronger at high redshift than at z — 0. This is because of 
the dramatic decrease of the non-linear clustering mass M* 
with redshift (see figure ^) . 

The time evolution of 6 var for halos (Af > 10 12 /t _1 M©) 
and galaxies (Mb — 51og/i < —19.5) is compared with the 
predictions of the models in figure Ell One can clearly see 
the slower evolution of gala xies versus halos in the simu- 
lations, 



discussed in Section 5.4. We show the MW model 
(eqn. [131) for the corresponding mass threshold. The MW 
model gives a good qualitative description of the evolution 
of halo biasing, as noted before. 

There is no mass threshold that gives a good description 
of the evolution of the galaxy biasing, because of the differing 
merger timescales of galaxies and halos, and because of the 
effects mentioned earlier concerning the changing mass-to- 
light connection as a function of redshift. 
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We also show the galaxy-c onserving (GC) mode l (|Fr; 
1996; Iegmark & Peebles 1998) described in Section 2.2.1. 
We start at high redshift [z = 3 or z — 1) with the values 
of 6 V ar and the linear correlation coefficient r as measured 
from the simulation galaxies or halos, and propagate the 
bias to z = using eqn. [ll]. The evolution predicted by the 
GC models is much weaker than that predicted by the MW 
models and observed in the simulations. This clearly sug- 
gests that merging, which is not included in the GC model, 
is an important process in the evolution of clustering over 
this redshift range, even in a low-f2 Universe. We might ex- 
pect that merging would be less important for galaxies than 
halos, and indeed the weak evolution according to the GC 
model is slightly closer to the biasing evolution for the galaxy 
population, but it still significantly underpredicts the rate 
of evolution. 



7 SUMMARY AND CONCLUSIONS 

We have performed a study of biasing of dark-matter halos 
and galaxies in cosmological simulations of ACDM (Q = 0.3) 
and rCDM (SI = 1), explicitly treating and quantifying the 
mean biasing as well as non-linearity and stochasticity in 
the biasing relation. We conclude that the CDM-based hi- 
erarchical structure formation scenario predicts that biasing 
has a moderate degree of nonlinearity and stochasticity, and 
it depends on mass, type, scale and especially redshift. 

The mean biasing function is always of the following 
characteristic shape. In the underdense regime (— 1 < 8 < 
0), b(8) vanishes near 8 £ — 1, then rises sharply towards 
<5 g ~ 8 = with a slope steeper than unity. In the overdense 
regime (8 > 0), the behavior is fairly linear, and the effective 
slope of b(S) is either larger or smaller than unity, represent- 
ing biasing or anti-biasing. The nonlinearity at 8/t _1 Mpc is 
at the level of a few to ~10 percent. The stochasticity is 
typically at the level of a few tens of percent. 

For halos, we find that the mean biasing increases grad- 
ually with mass, as does the stochasticity (to a large extent 
due to an increasing shot noise contribution) , while the non- 
linearity is nearly the same for all masses. At z = 0, the bias- 
ing relation for galaxies in our simulations is similar to that 
for > lO 12 /^ 1 Mq halos. This reflects the fact that our ha- 
los of 10 12 /i _1 Mq contain on average one bright galaxy per 
halo, but is still somewhat surprising given that the halos 
do contain different numbers of galaxies with various lumi- 
nosities. We find that the biasing is nearly independent of 
the luminosity of the galaxies selected, in contradiction with 
some observational results, but in keeping with our finding 
of weak mass dependence of halo biasing over the mass scales 
typical of galaxy-sized halos at z — 0. 

This apparent contradiction may be resolved by con- 
sidering dust extinction. We show that including the effects 
of dust effectively increases the mass-to-light ratio, so that 
L > L» galaxies are found in larger mass halos where the 
biasing dependence on mass is more pronounced. This leads 
to a luminosity dependent biasing that is qualitatively sim- 
ilar to that observed, although neither of the cosmological 
models considered here fully succeed in simultaneously re- 
producing the amplitude of the observed correlation func- 
tion and the detailed luminosity dependence of galaxy clus- 
tering. We suspect that a cosmology with an intermediate 



value of Qo ~ 0.5 might give better results. However, the 
observational results are still somewhat controversial, and 
the modelling of dust extinction is highly uncertain. Since 
this appears to have the potential to be a rather powerful 
constraint, the relative biasing of galaxies of different lumi- 
nosities should be measured more accurately using larger 
samples. Ideally, this should be investigated using a K-band 
selected sample which will decrease the uncertainties con- 
nected with dust extinction. 

We find that galaxy biasing depends fairly strongly on 
morphological type and color, with early type and red galax- 
ies being respectively about 1.3 to 1.8 times more biased 
than late type or blue galaxies. These results are in good 
agreement with current observational estimates. More de- 
tailed studies of relative biasing will be an important ap- 
plication of forthcoming large redshift surveys such as 2dF, 
SDSS, and 2MASS. 

Over the range of smoothing scales from 4 to 
16/i -1 Mpc, the mean biasing increases weakly with smooth- 
ing scale, similarly for halos and galaxies. 

The non-linearity and stochasticity decrease with in- 
creasing smoothing scale, so that the linear deterministic 
biasing approximation is approached for very large smooth- 
ing scales, as expected. On scales larger than a few Mpc, 
the clustering of bright galaxies in our simulations is fairly 
robust to differing assumptions about the details of galaxy 
formation. On smaller scales, halos become strongly anti- 
correlated because of exclusion effects. This implies that in 
order to reproduce the observed unbroken powerlaw corre- 
lation function of galaxies on these scales, there must be 
a significant number of galaxy pairs that cohabit the same 
halo. The biasing of our simulated galaxies on these scales 
is highly sensitive to the details of the astrophysical recipes, 
e.g., supernovae feedback, dust extinction, etc. 

The mean biasing of halos of a fixed mass and smooth- 
ing scale increases dramatically with redshift, by a factor of 
5-10 for rCDM and a factor of ~ 3 for ACDM from z = 
to z = 3. In both models, the nonlinearity of the biasing re- 
lation evolves with redshift, decreasing to a minimum value 
very close to 1.0 (pure linear bias) at z ~ 1 for rCDM and 
z ~ 3 for ACDM, then increasing again at higher redshift. 
The redshift dependence of biasing of galaxies with a fixed 
luminosity threshold differs significantly from that of halos 
with a fixed mass threshold. This is because, due to our 
standard star-formation recipe, galaxies of a given luminos- 
ity are hosted at higher redshift by halos of smaller masses 
and thus of weaker biasing. 

We have compared the results of the numerical simu- 
lations with analytic models for biasing. The MW model 
provides a good description of the changing shape of the 
biasing relation for halos as a function of mass threshold, 
smoothing scale, and more qualitatively, redshift. Even the 
simplified linear version of the MW model (Eqn. |l|) pro- 
vides a fairly accurate description of the mass dependence 
of mean halo biasing and the redshift evolution of halos of a 
fixed mass. The evolution of galaxy biasing may differ sub- 
stantially depending on the efficiency of star formation and 
how it va r ies with redshift. Galaxy-c onserving models (e.g. 
Fry 1996 ; Tegmark fc Peebles (1998 )) substantially under- 
predict the evolution rate of halo and galaxy biasing at all 
epochs, strongly suggesting that merging is a crucial ingre- 
dient in the evolution of biasing over this redshift range. 
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Our results are relevant to the attempts to measure the 
cosmological density parameter using galaxy densities, e.g., 
via redshift distortions or via a comparison to streaming 
velocities. Under the simplified assumption of linear and de- 
terministic biasing, these attempts lead to a measure of the 
quantity f3 = Q /b. Dekel & Lahav (1998) have proposed 
that deviations from the linear deterministic biasing ansatz 
may partially reconcile the differing values of (3 obtained by 
the different methods. Our results suggest that the expected 
levels of nonlinearity and stochasticity on the relevant scales, 
as predicted by simulations with realistic models of galaxy 
formation, would lead to modest differences in the various 
measures of j3, on the order of 20-30 percent. We suggest 
that the strong type dependence of biasing that is present 
both in our simulations and in observations may also con- 
tribute to some of the discrepancies in the results from dif- 
ferent surveys, which include different mixtures of galaxy 
types. This should be investigated further by applying these 
methods to detailed mock catalogs extracted from simula- 
tions similar to those used here. 

Our results also suggest several cautions that should be 
applied to the numerous recent attempts to draw conclu- 
sions from the redshift evolution of galaxy clustering. This 
may differ significantly from the redshift evolution of halo 
clustering because of differences in merging rates and a time 
dependence of the efficiency of star formation. Moreover, if 
different types of galaxy are selected at different redshifts, 
their biases may differ, giving a warped view of the actual 
redshift evolution. Finally, we find that different statistics 
for measuring the mean biasing may differ by as much as a 
factor of 2, and the disagreement of various statistics changes 
with epoch in accord with the changing importance of non- 
linearity, stochasticity, and scale dependence. These statis- 
tics are often treated as equivalent in the literature, an unfor- 
tunate outgrowth of the linear deterministic biasing ansatz. 

Our results are in general agreement both with other re- 
cent work using different methods and with what we know 
about the scale, type and redshift dependence of galaxy bi- 
asing from observations. Further studies using these tech- 
niques, along with new data from forthcoming large redshift 
surveys will no doubt lead to improved measurements of 
the cosmological parameters using galaxy-based methods, 
and in addition to a better understanding of the process of 
galaxy formation and evolution. 
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Figure 1. The relationship between galaxy luminosity and the mass of the host halo. The contours represent the joint probability that 
a galaxy has a given luminosity and dwells in a host halo with a given mass. The grey-scale indicates the log 10 of the probability, as 
shown by the scale bar on the figure. The right panels show the effects of applying a differential correction for dust extinction to the 
modelled luminosities, in which more luminous galaxies are more extinguished (see text). Top panels show the rCDM simulations and 
bottom panels show the ACDM simulations. The diagonal lines indicate constant mass-to-light ratio, with the constant dictated by the 
observed B-band Tully-Fisher relation. 
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Figure 3. The characteristic non-linear mass M, (see text) as 
a function of rcdshift for the tCDM and ACDM models. The 
horizontal dotted line indicates the mass of the smallest halos 
that we can reliably resolve in our simulations. 
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Figure 4. The probability of finding n galaxies per halo, with 
Mb — 51og/i < —19.5, as a function of halo mass, in the tCDM 
(top) and ACDM (bottom) simulations. The grey scale indicates 
the log 10 of the probability, according to the scale shown on the 
figure. The mean is shown by the cross symbols. 
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Figure 2. The mass to light ratio of halos in the simulations. The solid black line shows the B-band M/L, and error bars show the 
1-sigma scatter. The dotted line shows the B-band M/L after a correction for dust extinction has been applied (see text). The grey line 
shows the I-band M/L (without dust). 
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Figure 5. The joint distribution of the overdensity fields of halos (M > 10 1 h^ 1 Mq) and mass, both smoothed with a top-hat window 
of radius R B = 8h~ 1 Mpc, for rCDM (left) and ACDM (right), for several redshifts. The contours represent approximately the 50 (dark 
grey), 80 (medium grey), and 98 (light grey) percentiles. The white lines show the mean conditional biasing function (5^|<5) and the la 
scatter about it in equal bins of log(l + <5). The black long-dashed line shows the linear biasing approximation with b = b. 
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Figure 7. Mass dependence of the mean and scatter of the condi- 
tional biasing relation for halos, for rCDM at z = 0. The bottom 
panel shows the mean conditional biasing relation for halos in 
our simulations, selected with different mass thresholds. The top 
panel shows the variance of the conditional biasing relation, minus 
the mean expected variance due to shot noise. Bold lines indicate 
the simulation results and light lines indicate the "conditional 
Poisson" model discussed in the text. 



Figure 8. Same as figure for galaxies selected with different 
B-band absolute magnitude thresholds. 
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Figure 10. Biasing of halos as measured by the two-point auto- 
correlation functions for different halo-mass thresholds at z = 
in the rCDM simulations. The bold solid line refers to the mass 
and the broken lines to the halos. The fit to the observed galaxy 
correlations in the APM survey is shown for reference (light solid 
line). Results for the ACDM simulations are qualitatively similar. 



Figure 9. Biasing characteristics as a function of halo mass 
threshold and galaxy magnitude limit. The simulations are an- 
alyzed at z = 0, with T8 smoothing. Top: mean biasing 6, 6 var , 
and fej. Middle: nonlinearity b/b. Bottom: stochasticity cr^/b. 
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Figure 11. The auto-correlation function for galaxies (light broken lines) with varying absolute magnitude limits at z = in the rCDM 
simulation (top) and ACDM simulations (bottom). The right panels include a model for dust extinction (see text). Symbols sho w the 
correlation length j-q for galaxies with varying absolute magnitude thresholds, obtained by Willmer, Da Costa, & Pellegrini (1998) from 
the SSRS2 redshift survey. 
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Figure 12. The joint distribution of the overdensity fields of type-selected galaxies and all galaxies, both with Mg — 51og/i < —18.4 
and smoothed with a top-hat window of radius R s = 8/i _1 Mpc, for the rCDM simulations. The top panels represent galaxies selected 
by morphology ("early" corresponds roughly to E-S0, "late" to S-Irr Hubble types), and the bottom panels represent galaxies selected 
according to B-V colour. The contours represent approximately the 50 (dark grey), 80 (medium grey), and 98 (light grey) percentiles. 
The white lines show the mean conditional biasing function ((Jeariy/iatel^g) an d tne la scatter about it in equal bins of log<5 g . The black 
long-dashed line shows the linear biasing approximation with b = b rc i . 
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Figure 13. The distribution of host halo masses for galaxies selected according to morphological type (left) and colour (right) and with 
Mg — 51og/i < —18.4. Solid lines show the distribution of host halo masses for all galaxies above the magnitude limit, dotted lines 
show the distribution of masses for late/blue galaxies, and dashed lines show the distribution for early/red galaxies. Early-type and red 
galaxies dwell on average in more massive halos, and hence are biased compared to the overall galaxy population. 
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Figure 14. Scale dependence of the mean and scatter of the Figure 15. Same as figure |14|, for galaxies with Mb — 51ogh < 

conditional biasing relation for halos (M > 10 12 /i — 1 Mq), for —19.5. 
tCDM at z = 0. The quantities shown are as in figure [7| 
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Figure 16. Biasing characteristics as a function of smoothing 
scale for halos of M > I0 12 h~ 1 Mq or galaxies with Mg — 
51og/i < —19.5 at z = 0. Details are as in figure Bl 



Figure 17. Biasing characteristics as a function of redshift. The 
two cosmological models are analyzed with T8 smoothing, for 
halos of M > lO 12 /!" 1 Mq or galaxies with M B -5 log h < -19.5. 
Details are as in figure M. 
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Figure 18. The auto-correlation function for halos (M > I0 12 h 1 Mq; left panels, light lines) or galaxies (Mg — 51og/i < —19.5; right 
panels, light lines) and for the dark matter (bold lines) at different rcdshifts. 
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Figure 20. Mass dependence of halo biasing (fe V ar) at z = 0, 
z = 1, and 2 = 3. Solid lines are the predictions of the MW model. 
Dotted lines arc results from the tCDM simulations without cor- 
rection for shot-noise. Symbols are results from the simulations 
with a standard Poisson correction for shot-noise. Circles and the 
lowest set of lines arc for 2 = 0, diamonds and the middle set of 
lines for 2 = 1, and triangles and the highest lines are for z = 3. 



Figure 19. Comparison of the predictions of the MW model 
(light lines) with the simulation results (bold lines) for the mean 
biasing relation. The curves in each panel have been offset by 
0.5 dex for clarity. The top panel shows the results for different 
mass thresholds (loglM/^- 1 M )] = 11.4,12.0,12.4,12.7 from 
left to right), the middle panel shows different smoothing scales 

(R s = 4,8, 12, Wh,- 1 Mpc from left to right) and the bottom © 0000 RAS, MNRAS 000, 000-000 

panel shows different rcdshifts (2 = 0, 1, 2, 3 from left to right). 
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Figure 21. Redshift evolution of biasing (6 V ar) of galaxies 
(M B - 5 log h < -19.5; stars) and halos (M > lO 12 /*^ 1 M ; 
dots) in the simulations, compared with analytic models. The 
light solid lines show the predictions of the MW model, with the 
same mass threshold as the simulation halos. The dotted lines 
show the predictions of the galaxy-conserving model, for initial 
bias values typical of halos or galaxies at z = 3 or galaxies at 
z = 1 in our simulations (see text for details). 
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